Chirality of High Latitude Filaments over Solar Cycle 23 
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ABSTRACT 

A non-potential quasi-static evolution model coupling the Sun's photospheric 
and coronal magnetic fields is applied to the problem of filament chirality at high 
latitudes. For the first time, we run a continuous 15 year simulation, using bipolar 
active regions determined from US National Solar Observatory, Kitt Peak mag- 
netograms between 1996 and 2011. Using this simulation, we are able to address 
the outstanding question of whether magnetic helicity transport from active lat- 
itudes can overcome the effect of differential rotation at higher latitudes. Acting 
alone, differential rotation would produce high latitude filaments with opposite 
chirality to the majority type in each hemisphere. We find that differential rota- 
tion can indeed lead to opposite chirality at high latitudes, but only for around 5 
years of the solar cycle following the polar field reversal. At other times, including 
the rising phase, transport of magnetic helicity from lower latitudes overcomes 
the effect of in situ differential rotation, producing the majority chirality even 
on the polar crowns at polar field reversal. These simulation predictions will 
allow for future testing of the non-potential coronal model. The results indicate 
the importance of long-term memory and helicity transport from active latitudes 
when modeling the structure and topology of the coronal magnetic field at higher 
latitudes. 



Subject headings: Sun: corona — Sun: evolution — Sun: filaments, prominences — 
Sun: magnetic topology — Sun: surface magnetism 
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Introduction 



Observations such as X-ray sigmoids, flares, or eruptions strongly indicate the presence 
of substantial electric currents in the low solar corona (below, say, 2.5it^), which are ab- 



sent by definition in traditional potential field extrapolations (lAltschuler fc Newkirkl 11969 



Schatten et al.lll969[ ). To account for the development and effects of these electric currents, 
we have introduced a non- potential (hereafter NP) model coupling the Sun's pho tospheric 
and coronal magnetic fields ( Ivan Ballegooijen et al.ll2000l ; lYeates et al.ll2008al . l2010l ). The NP 
model circumvents the problem of non-uniqueness which plagues the e xtrapolation of non- 
linear force-free equilibria from observed photospheric magnetograms ( iDerosa et al.l 120091 ) 
by building up currents in a time-dependent manner. The coronal magnetic field relaxes 
continually (toward a nonlinear force- free state) in response t o shearing by large-scale pho- 
tospheric motions, and new flux emergence. Previous testing (lYeates et al.ll2008al |2010| ) has 
shown this model to better approximate the real physics determining the large-scale coronal 
magnetic field. 



To test the global NP model, lYeates et al.l (j2008al ) compared the simulated magnetic 
field direction with the chirality of quiescent solar filaments observed in Ha over a 6-month 
period, finding excellent agreement. Solar filaments are stable structures of (relatively) cool, 



dense plasma suspended well above the chromosphere ( iLabrosse et al.l l2010t iMackay et al. 



20101 ) . and overlying Polarity Inversion Lines (PILs) in the photospheric magnetic field. They 
are seen in both Ha and Extreme Ultraviolet, either against the solar disk or as prominences 
above the limb. The large-scale magnetic field in filaments is believed to be orientated along 
the filament axis. Filaments may then be classified as having either dextral or sinistral 
chirality, according to the direction of this axial field relative to the underlying photospheric 



polari ties (Figure [JJ. Although a property of the filament magnetic field, iMartin et al. 



( 119941 ) showed that the chirality may often be determined from the orientation of barbs in 
Ha images, given only the photospheric polarities on either side of the filament channel. 
Observations reveal a stati stical preference for dextral chirality in the northern hemisphere 



and sinistral in the south (Rust 



19671 iLerov et al.l 11983 : Martin et all 11994 : IPevtsov et al. 



20031 ; iBernasconi et al.ll2005l ). as indicated in Figure (Ha) 



For filaments equatorward of about 60° latitude, we were able to use the NP model to 
explain both the observed pattern and how exceptions can occur. The chirality results from a 
number of c ompeting effects, depending on the fi lament location ( IMackay fc van Ballegooijen 
200 ll . hood : lYeatesll2009l ; lYeates k Mackavll2009h . PILs with in individual active r egions have 
chirality consistent with the helicity sign in the active region ( iRust fc Martinlll994l ). But PILs 
between or outside of active regions can be strongly affected by interactions between multiple 
regions — this explains many exceptions to the majority hemispheric pattern (supported by 
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the observations of lWang fc Muglachll2007l ; lRomano et al.ll2009l ). Also, over ti me, differentia l 
rotation can strengthen, weaken, or even reverse the chirality at a given PIL (jYeatesll2009l ). 



At higher latitudes, projection effects make it difficult to determ ine filament chiralit y 
through Ha barbs (for example, polar crown filaments were omitted by lPevtsov et al.ll2003l ). 
However, several sets of linear polarization measurements in limb prominences indic ate the 
same hemispheric chirality patter n, both at solar m inimum (jRustlll967l ; iLeroy Ill978l ) and in 
the rising phase of the solar cycle (ILeroy et al.lll983l ). although (to our knowledge) no studies 
during the declining phase have been published. This finding that the same hemispheric 
pattern holds at higher latitudes was not reproduced by our previous 6-month simulations 
for 1999 (in the rising phase of Cycle 23), where the polar crown chirality was found to be 
opposite to that of lower latitudes. 

We attributed this result to the dominance of differential rotation on the East- West 
oriented polar crown PILs. As shown in Figure Hf b), differential rotation alone l eads to the 
wrong chirality overlying East- West oriented PILs (ILeroy Hl978l ; IZirker et al.lll997l ). However, 
the chirality of high-latitude PILs in the NP model is a balance between differential rotation 
and the poleward transport of magnetic helicity from lower latitudes. (We are assuming that 
there is no net injection of helicity from the unresolved small-scale motions.) The possibility 
therefore remains that there was simply insufficient time in the 6-month simulations for 
poleward transport of helicity. In this Letter we extend the simulation to 15 years' duration, 
simulating both the rising and decaying phases of the solar cycle, to explore the predicted 
chirality at high latitudes. 



2. Simulation of Cycle 23 



We simulate the evolution of the large-scale photospheric and coronal magnetic fields 
between 1996 and 2011, using the non-potential (NP) model , which couples surface flux 



t r ansp ort to magnet o-fric t ional relaxation in the corona above (Ivan Ballegooijen et al.ll2000 



Yeates et al.ll2008al . |2010| ). The large-scale mean magnetic field, B 
is evolved via the induction equation 



V x A , in the corona 



dt 



v x B - E 



where we neglect Ohmic diffusion and the mea n electromotive force En descr i bes th e effect of 
unresolved small-scale fluctuations. Following Ivan Ballegooijen Sz Cranmerl (120081 ) we apply 
a hyper diffusion 



En = -^V- 



' r] 4 Bl V« 



Bl 



(2) 
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(a) Observed chiralities 




(b) Expected chiralities from 
differential rotation alone 



sinistral 




dextral 



Fig. 1. — Filament magnetic fields are assigned dextral or sinistral chirality according to 
the angle between their horizontal field and the unde rlying PIL. Observations indicate the 



global chirality pattern in (a) (e.g.. iLeroy et al.lll983l ). while di fferential acting a lone on an 



unsheared magnetic arcade would lead to the pattern in (b) ( jZirker et al.lll997l ). See the 
electronic edition of the Journal for a color version of this figure. 
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where 



B ■ V x B 



a 



o 



Bl 



(3) 



and we take 774 = 10 11 km 4 s~ 1 . This form of hyperdiffusion preserves magnetic helicity 
density Ao • Bo in the volume and describes the tendency of the m agnetic field to relax to 
a state of constant «o ( lBoozerlll986t iBhattacharjee fc Hameirilll986l ). although such a state 
is never reached in the global simulation. We have repeated the simulation with a constant 
diffusivity (Eo = — 77Jo) with no cha nge to the chiralit y results. The velocity is determined 
by the magneto- frictional technique flYang et al.lll986l ) as 

1 jo x B , v , s 

v = + v out {r)r, (4) 

where the first term enforces relaxation toward a force-free equilibrium, and the second term 
is a radial outflow imposed only near the outer boundary (r = 2.5R ( r ) ) to represent the effec t 
of the solar wind radially distending magnetic field lines ( iMackay fc van Ballegooijenll2006l ). 



To run the simulation continuously over 15 years, we use a non-uniform spherical grid 
whose resolution reduces from 192 grid points in longitude at the equator to only 12 at 
the polar grid boundaries (±89.5° latitude). The simulation is initialised with a potential 
field extrapolation based on a synoptic magnetogram for Carrington Rotation CR1910 from 
US National Solar Observatory, Kitt Peak (NSO/KP), corrected for differential rotation to 
represent 1996 May 15. The coronal magnetic field is then evolved continuously, driven 
by the emergence of new active regions and shearing by large-scale photospheric motions. 
The latter are im posed on the lower boundary r = R & through the standard surface flux 



transport model (ILeighton 



D = 450km V\ the 



1964 



Snodgrass 



Sheelevi 120051 ) . Here we use a supergranular diffusivity 



1982 ) differential rotation profile, and the meridional flow 
profile of ISchiissler fc Baumannl (120061 ) with peak speed 11 ms _1 . In order for the simulation 



to reproduce the low polar field stre ngths observed after r eversal in Cycle 23, we reduce 
the active region tilt angles by 20% (ICameron et al.ll2010l ; I Jiang et al.l l201ll ). and include 
an additional decay term for the photosph eric magnetic field with a timescale of 10 years 
(iSchrijver et alihood : ISchrijver fc LkJhoQsl ). 



Active regions a re inserted in the form of idealised three-dimensional magnetic bipoles 
(jYeates et al.l l2008al ). whose properties (location, size, magnetic flux and tilt angle) are 
determined from NSO/KP synoptic magnetograms. Until 2003 these were taken with the 
older vacuum telescope, and from 2003 onwards with SOLIS (Synoptic Optical Long-Term 
Investigations of the Sun). For the present simulations we do not insert any bipoles to replace 
those missed during the three data gaps (CR2015-16, CR2040-41, CR2091). In total, between 
CR1911 (1996 June) and CR2110 (2011 May), we insert 1838 bipoles. The regions vary in 
flux content from 2 x 10 20 Mx to 5.3 x 10 22 Mx. 
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As stated above, we reduce the tilt angles by 20% to help reproduce the low polar fields 
observed after reversal in Cycle 23. A linear fit to the reduced tilt angles gives sin(tilt) = 
0.38 sin(Ao)— 0.005, wh ich reduces the slope of t he relationship compared to the magnetogram 
observations of either IWang fc Sheeleyl (119891 ) (who used NSO/KP data for Cycle 21) or 
Stenflo &: Kosovichevi ( 20121 ) (who used SOH O/MDI magnetog r ams) . However, it is still 
greater than found in the white light study of iDasi-Espuig et al.l ( 120101 ) for Cycles 15 to 21. 
In principle, the same resulting polar fields could be obtained wit h the original bip ole tilt 
angles by instead invoking a more complex meridional flow profile (1 Jiang et al.ll2010l ). 



The three-dimensional bip oles have a twist pa rameter /3 describing emergence of mag- 
netic helicity within the bipole (jYeates et al.ll2008al ). Since the vector magnetic observations 
required to determine f3 are not available for ev ery active region in Cycle 23, we select (5 statis- 
tically for each new bipole ( jYeates et al.ll2010l ). Although we have s hown that this emergin g 
helicity affects the chirality of ~ 32% of filaments at lower latitudes (jYeates &: Mackayll2009l ). 
these are mainly those located within recently-emerged active regions. After only 30-50 days, 
differential rotatio n tends to rev erse the chirality of PILs inside any bipoles inserted with the 
minority sign of /3 (jYeatesll2009l ). Thus the choice of /3 distribution has little consequence for 
the higher latitudes considered here. We have verified that the same high latitude chirality 
pattern is recovered if ft — for all bipoles. 



3. Results 

Figure [2] shows that the photospheric radial magnetic field B Qr in the observed magne- 
tograms (Figure^) is generally well reproduced by the simulation (Figure [2b). 

To analyse the chirality, snapshots of the global three-dimensional magnetic field were 
stored every 14 days during the simulation. For each snapshot, PILs were identified automat- 
ically on the photosphere r = R Q , along with the overlying chirality of the coronal magnetic 
field at height r = r\. We take T\ = 1.033i? Q to approximate the height of filaments. Now 
denote the horizontal field at r — T\ by B^, and the radial field in the photosphere by B 0r . 
We characterise the chirality at height r = n by the "skew" 

e r • VB 0r x Bf 

sm 7 = — i n — ri — ' w 

so that sin 7 < corresponds to sinistral skew and sin 7 > to dextral skew. Two example 
snapshots from the simulation, on 1999 June 24 and 2002 June 20, are shown in Figure |3j 
The color shading along photospheric PILs indicates the skew calculated above them. 



Figure Ufa) shows the longitude- averaged skew angle along PILs in 95 latitude bins, as 
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Year 

Fig. 2. — Butterfly diagrams showing longitude-averaged radial magnetic field Bq t over time, 
in (a) NSO/KP synoptic magnetograms and (b) the simulation. Three data gaps are evident 
in (a). 
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Fig. 3. — Snapshots of the Cycle 23 non-potential simulation (a) before and (b) after polar 
field reversal. Left images show photospheric PILs colored according to skew sin 7 at r = 
l.O33i? (in Mollweide equal-area projection), while right images show selected coronal field 
lines traced from this height. Dextral skew corresponds to sin 7 > (red) and sinistral 
to sin 7 < (blue). In both cases, grey shading shows B 0r on the photosphere R & (white 
positive, black negative; saturation ±30 G). The left images are also available as an mpeg 
animation in the electronic edition of the Journal. 
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a function of time through the simulation. Several features are apparent: 



1. 



In and around active regions the majority pattern of chirality predominates, i.e., dextral 
in the northern hemisphere and sinistral in the southern hemisphere, although this is 
an overall mean pattern with significant local fluctuations in both magnitude and sign 
(as seen in Figure El) . The o rigin of t his pattern, in cluding the exceptions, is explained 
bv lYeates fc Mackavi f|2009h (see also lYeates^2009h . 



2. During the period of few active regions from 2007 to 2010, there is more mixed chirality 
at lower latitudes. 

3. Until 1998, and during the declining phase from 2001 to 2006, there is a tendency 
for minority chirality on the high-latitude PILs (sinistral in the north, dextral in the 
south). 

4. During the "rush-to-the-poles" between 1998 and mid-1999, the polar crowns exhibit 
the majority chirality pattern. 

5. From 2006 onward, the majority chirality dominates at high latitudes once more, con- 
tinuing into the start of Cycle 24. 



For co mparison, HKc) show s the longitude-averaged current helicity density «o (Equation 
El see also lYeates et al.ll2008bl ). at the same height r = 1.033f? Q . As expected, the sign of 
«o correlates approximately with that of the skew angle, where positive ao is associated 
with negative skew angle, and vice versa. (The concentrations of ao at the poles are due to 
twisting up of the vertical field there by the differential rotation.) 



Discussion 



The chi rality pattern at hig h latitudes in Figure H(b) differs from the earlier 6-month 
simulations (lYeates et al.l l2008a| ) . In those simulations, high-latitude filaments remained 
predominantly sinistral in the northern hemisphere and dextral in the southern hemisphere: 
opposite to the majority pattern. In the new 15-year simulation, the majority pattern is 
found in two independent epochs: once during the rising phase up to polar reversal (including 
the period of the previous simulations), and again from the middle of the declining phase of 
Cycle 23. We argue that: 



1. Transport of helicity from active latitudes is responsible for overcoming the effect of 
differential rotation on the polar crown over much (but not all) of the cycle. 
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Year 

Fig. 4. — Butterfly diagrams showing longitude-averages of (a) skew sin 7, and (b) current 
helicity density cxq, during the Cycle 23 non-potential simulation. Both are measured at 
height r = l.O33i? . See the electronic edition of the Journal for a color version of this 
figure. 
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2. The initial period of opposite chirality (before 1999) is an artefact of the initial condi- 
tion of the simulation, indicating the importance of magnetic memory in the corona. 
Note that such a period of opposite chirality does not recur at the start of Cycle 24 
(from 2009). 

An example where differential rotation clearly wins over helicity transport is shown in 
Figure El^b), during the declining phase of the solar cycle. Here, minority dextral chirality is 
found on the poleward arm of a switchback in the southern hemisphere. The field lines (Fig- 
ure [31 bottom right) show a sheared arcade, long disconnected from active region flux. The 
tendency for minority chirality on such switchbacks in the declini ng but not in the rising phase 



match es the predictions made from single-bipole simulations by iMackay fc van Ballegooijen 



( 1200 if ). 



By constrast, the north polar crown during the rising phase prior to polar reversal (Fig- 
ure 13^) exhibits a strong flux rope with majority chirality, with field lines directly connecting 
to lower latitudes. To further demonstrate that the polar crown chirality is controlled by 
a balance between differential rotation and helicity transport, we repeated the simulation 
with the latitudinal differential gradient reduced to half of its original slope. This leads to 
less minority chirality at high latitudes (Figure |5]), consistent with the greater dominance 
of helicity transport. Even in the original simulation, majority chirality takes over at all 
latitudes later in the cycle (from 2006 onward), because the highest-latitude PILs at this 
time are nearer to the equator than those earlier in the cycle, so that helicity from lower 
latitudes can more easily dominate. 

The initial period of minority chirality at high latitudes does not recur at the start of 
Cycle 24 in the simulation. This suggests it may be an artefact of the initial potential field. 
Majority chirality on the polar crowns during the minimum and rising phase of the solar 
cycle w ould be requ i red in order to be consistent with the reported polarization measure- 



ments (ILeroy I Il978l ; iLeroy et al.lll983l ). Indeed, when we simulate two consecutive cycles 
with emerging active regions chosen at random from statistical distributions (not illustrated 
here), we find this minority chirality only at the start of the simulation, not at the start 
of subsequent cycles. Because the high-latitude field depends on transport of helicity from 
lower latitudes, differential rotation dominates early in the simulation, producing minority 
chirality on the polar crowns. In effect, the high latitude field has a finite memory of around 
12 months for the pre-existing magnetic topo logy at the start of the simulation. This ex- 



plains why the previous 6-month simulations ( lYeates et al.ll2008al ) found minority chirality 
on the polar crowns in 1999. 

In conclusion, the NP model makes robust predictions for the chirality of high-latitude 
filaments: majority chirality over much of the cycle (including during polar field reversal), 
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but more mixed chirality between the polar reversal and middle of the declining phase. The 
first prediction matches the available polarization measurements of prominences from earlier 
solar cycles, but further measurements are needed to confirm the second prediction. Mixed 
chirality is also predicted at lower latitudes during particularly quiet periods such as the 
recent Cycle 23 minimum. There are some indications that the hemisphe ric pattern of actiy e 
region helicity was also weaker during the d eclining; phase of Cy cle 23 (ITiwari et al.ll2009l ). 
although this remains statistically uncertain ( iPevtsov et al.ll2008l ). However, the multitude of 
factors affecting filament chirality mean that we do not expect to see the same time variation s 
in the filament and active region hemispheric patterns (lYeatesll2009l ; lYeates &: Mackayll2009l ). 
Testing our predictions against prominence measurements will best validate the NP model 
for the global coronal magnetic field. Of course, our model at present considers only the 
large-scale magnetic field within which filaments form. The smaller scale magnetic field 
structure inside filaments /prom inences themselves is inferred to be rather complex and is 
the subject of ongoing research (IMerenda et al.ll2006l ; Ivan Ballegooijen &: Cranmerll2010l ). 
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Fig. 5. — Longitude-averaged skew sin 7 at height r = 1.033i? Q for the simulation with 
reduced differential rotation, showing reduced areas of minority chirality. (Note that the 
longitu de-averaged Bn r in the photosphere is not affected by changing the differential ro- 
tation; 
figure. 



Leightonlll964l ). See the electronic edition of the Journal for a color version of this 



